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Abstract. The phonon dynamics of the low-temperature superconductor Sr2Ru04 is 
calculated quantitatively in linear response theory and compared with the structurally 
isomorphic high-temperature superconductor La2Cu04. Our calculation corrects for 
a typical deficit of LDA-based calculations which always predict a too large electronic 
fcz-dispersion insufficient to describe the c-axis response in the real materials. With 
a more realistic computation of the electronic band structure the frequency and 
wavevector dependent irreducible polarization part of the density response function 
is determined and used for adiabatic and nonadiabatic phonon calculations. Our 
analysis for Sr2Ru04 reveals important differences from the lattice dynamics of p- 
and n-doped cuprates. Consistent with experimental evidence from inelastic neutron 
scattering the anomalous doping related softening of the strongly coupling high- 
frequency oxygen bond-stretching modes (OBSM) which is generic for the cuprate 
superconductors is largely suppressed or completely absent, respectively, depending 
on the actual value of the on-site Coulomb repulsion of the Ru4d orbitals. Also the 
presence of a characteristic Ai-mode with a very steep dispersion coupling strongly 
with the electrons is missing in Sr2Ru04. Moreover, we evaluate the possibility of a 
phonon-plasmon scenario for Sr2Ru04 which has been shown recently to be realistic 
for La2Cu04. In contrast to La2Cu04 in Sr2Ru04 the very low lying plasmons are 
overdamped along the c-axis. 
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1. Introduction 

The discovery of superconductivity in Sr2Ru04 by Maeno and coworkers pLj has attracted 
widespread attention partially because of the structural similarity with the cuprate- 
based high-temperature superconductors (HTSC's). While the mother compounds of 
the HTSC's are charge transfer insulators, and they usually need to be doped to become 
metalhc and show superconductivity, in Sr2Ru04 low-temperature superconductivity 
at Tc = 1.5 K condenses from a metallic state that is a strongly two-dimensional 
Fermi liquid below about 20 K. The Fermi surface consists of three weakly corrugated 
cylindrical sheets a being hole-like and /3 and 7 which are electon-like p]. This is in 
contrast to La2Cu04 where only one Fermi sheet exists. The mechanism of conduction 
in Sr2Ru04 at higher temperatures is an interesting issue because the c-axis resistivity 
shows a broad maximum at around 130 K and for increasing temperatures the resistivity 
starts to decrease. In [3] such a metallic to nonmetallic crossover in the c-axis resistivity 
of a highly anisotropic metal like Sr2Ru04 with no corresponding feature in the afe-plane 
properties has been related to the strong coupling between the electrons and a bosonic 
mode propagating and polarized in c-direction, like the axial oxygen breathing mode 
Of, at the Z point of the Brillouin zone (BZ) to be discussed in section [3l 

Superconductivity in Sr2Ru04 is thought to be of unconventional character and 
thus is intensively discussed in the literature. The current state of the experiments 
points to triplet superconductivity in Sr2Ru04 in contrast to the d-wave singlet pairing 
in the HTSC's. Most discussion of the superconducting mechanism for Sr2Ru04 is 
focussed on magnetic fluctuations and Coulomb repulsion, for a review see [4J. However, 
a contribution of electron-phonon interaction (EPI) to pairing may not be ruled out 
because the superconducting transition exhibits a clear isotope effect |5j. In case of 
the HTSC's there is increasing evidence that EPI is strong and phonons might play an 
important role for the electron dynamics, see [6], [7] and references therein. 

Thus, a detailed theoretical study of phonon dynamics in Sr2Ru04 and a comparison 
with the situation in the HTSC's is desirable. The low-temperature superconductor 
Sr2Ru04 crystallizes in the K2NiF4 structure (I4/mmm) isostructural to the HTSC 
La2Cu04. An analysis and comparison of the lattice dynamics in both materials is the 
main intention of the present work. For this purpose we compute the phonon dynamics 
for Sr2Ru04 and oppose it with our earlier calculations for La2Cu04. 

The most characteristic feature of phonon dynamics in the HTSC's is the strong 
softening upon doping of the high-frequency oxygen bond-stretching modes (OBSM) 
which seems to be generic for p- as well as n-doped cuprates, see [El [3, [HI [9] . The strong 
frequency renormalization and a related increase of the linewidths observed for these 
anomalous modes point to a strong coupling of the phonons to the charge dynamics. A 
discussion and theoretical investigations of these phonon anomalies are presented e.g. 
in [HI [71 [HI [TOj. In agreement with inelastic neutron scattering [11] our calculations 
demonstrate that the typical OBSM phonon anomalies of the cuprates are not present 
in Sr2Ru04. Moreover the characteristic Ai branch related to the strongly coupling Of 
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mode in La2Cu04 is missing. 

In order to achieve reliable results for the Ai phonon modes propagating and 
polarized along the c-axis a careful calculation of the proper polarization part of the 
density response function is crucial. This has already been shown for La2Cu04 in 
[12] where we have proven that the Ai modes are highly sensitive with respect to 
the charge response perpendicular to the CuO plane. The large anisotropy of the 
electronic structure of the cuprates is as a rule considerably underestimated in DFT- 
LDA calculations and as a consequence the Ai modes are not well described. So we 
have modified a LDA-based calculation for La2Cu04 to account for the much weaker 
fc^-dispersion of the electronic bandstructure in the real material. 

We have optimized the interlayer coupling in such a way that the significant features 
of the sensitive Ai phonons are well described. The same route is followed in the present 
work for Sr2Ru04 to get a realistic representation of the Ai modes and in particular of 
of which cannot be described in the ionic shell-model taking additionally homogeneous 
electron gas screening into account [TTj . 

Finally, we investigate the question if there is room for a phonon-plasmon scenario 
around the c-axis for Sr2Ru04 which has been shown to be a realistic option in La2Cu04 
|12] . Due to the much weaker electronic fc^-dispersion obtained in our computations as 
compared to La2Cu04, i.e. a much stronger anisotropy of Sr2Ru04, the calculated free- 
plasmon frequencies along the c-axis are about a factor of eight smaller in the coUisionless 
regime. Damping, generated by interactions between the electrons as well as interband 
transitions which are on a much lower energy scale as for La2Cu04 very likely leads to 
overdamping of the plasmons propagating strictly along the c-axis in Sr2Ru04. Thus, 
different from La2Cu04 coupled c-axis phonon-plasmon modes should not be well defined 
collective excitations at least strictly along the A ~ (0, 0, 1) direction and accordingly 
do not enter the list of possible players for pairing in contrast to La2Cu04. 

The article is organized as follows. In section [2] the theory necessary to understand 
the calculated results is shortly reviewed. Section [3] presents the calculations. The 
modification of a LDA-based electronic bandstructure providing the basis of the single 
particle content of the irreducible polarization part of the density response function is 
developed. Furthermore, the phonon dispersion is calculated in adiabatic approximation 
for Sr2Ru04 and a comparative discussion with the situation in La2Cu04 is given. 
Finally, in a nonadiabatic calculation the possibility of a phonon-plasmon scenario is 
examined. A summary of the paper is given in section H] and the conclusions are drawn. 

2. Theory and modeling 

In the following a brief survey of the theory and modeling is presented. A detailed 
description can be found in [10] and in particular in [ISj where the calculation of the 
coupling parameters of the theory is presented. 

The local part of the electronic charge response and the EPI is approximated in the 
spirit of the quasi-ion approach [3, [H] by an ab initio rigid ion model (RIM) taking into 
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account covalent ion softening in terms of (static) effective ionic charges calculated from 
a tight-binding analysis. The tight-binding analysis supplies these charges as extracted 
from the orbital occupation numbers of the fi (tight-binding) orbital in question: 

Q, = ^Y.\^.n{k)\'. (1) 
nk 

C^„(k) stands for the /i-component of the eigenvector of band n at the wavevector k in 
the first BZ; the summation in ([1]) runs over all occupied states and gives the number 
of the elementary cells in the (periodic) crystal. 

In addition, scaling of the short-ranged part of certain pair potentials between the 
ions is performed to simulate further covalence effects in the calculation in such a way 
that the energy-minimized structure is as close as possible to the experimental one 
|15j . Structure optimization and energy minimization is very important for a reliable 
calculation of the phonon dynamics through the dynamical matrix. Taking just the 
experimental structure data as is done in many cases in the literature may lead to 
uncontrolled errors in the phonon calculations. 

The RIM with the corrections just mentioned then serves as an unbiased reference 
system for the description of the HTSC's and can be considered as a first approximation 
for the insulating state of these compounds. Starting with such an unprejudiced rigid 
reference system non-rigid electronic polarization processes are introduced in form 
of more or less localized electronic charge-fluctuations (CF's) at the outer shells of 
the ions. Especially in the metallic state of the HTSC's the latter dominate the 
nonlocal contribution of the electronic density response and the EPI and are particularly 
important in the CuO planes and the RuO plane. In addition, anisotropic dipole- 
fiuctuations (DF's) are admitted in our approach [13l[l6], which prove to be specifically 
of interest for the ions in the ionic layers mediating the dielectric coupling and for the 
polar modes. Thus, the basic variable of our model is the ionic density which is given 
in the perturbed state by 

Pa(r, Qa, Pa) = Plir) + Q^PTir) + Pa " fp°(r). (2) 

A 

p° is the density of the unperturbed ion, as used in the RIM, localized at the sublattice 
a of the crystal and moving rigidly with the latter under displacement. The Qx and 
describe the amplitudes and the form-factors of the CF's and the last term in equation 
([2]) represents the dipolar deformation of an ion a with amplitude (dipole moment) 
and a radial density distribution p^. r denotes the unit vector in the direction of r. The 
p^^ are approximated by a spherical average of the orbital densities of the ionic shells 
calculated in LDA taking self-interaction effects (SIC) into account. The dipole density 
is obtained from a modified Sternheimer method in the framework of LDA-SIC |13j . 
All SIC-calculations are performed for the average spherical shell in the orbital-averaged 
form according to Perdew and Zunger [17]. For the correlation part of the energy per 
electron the parametrization given in [T7j has been used. 

The total energy of the crystal is obtained by assuming that the density can be 
approximated by a superposition of overlapping densities p^. The p° in equation ([2]) 
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are also calculated within LDA-SIC taking environment effects, via a Watson sphere 
potential and the calculated static effective charges of the ions into account. The 
Watson sphere method is only used for the oxygen ions and the depth of the Watson 
sphere potential is set as the Madelung potential at the corresponding site. Such an 
approximation holds well in the HTSC's [T5[[18]. Finally, applying the pair-potential 
approximation we get for the total energy: 

EiR,0 = J2E^o.iC) + l E -^-a, (R' - R^, C) • (3) 

a,a {a,a)^(h,l3) 

The energy E depends on both the configuration of the ions {R} and the electronic 
(charge) degrees of freedom (EDF) {(} of the charge density, i.e. {Q\} and {pa} in 
equation ([2]). are the energies of the single ions, a, b denote the elementary cells and 
a, P the corresponding sublattices. The second term in equation ([3]) is the interaction 
energy of the system, expressed in terms of anisotropic pair- interactions $q,/3. Both E^ 
and in general depend upon ( via pa in equation ([2]). 

The pair potentials in equation ([3]) can be seperated into long-ranged Coulomb 
contributions and short-ranged terms, for details see e.g. Ref |13j . 

The dynamical matrix in harmonic approximation can be derived as 



RIM 



-^i:i«r(q)l'r'(q)]..i^r(q). (4) 

The first term on the right hand side denotes the contribution from the RIM. M^, 
are the masses of the ions and q is a wave vector from the first BZ. 

The quantities B(q) and C(q) in equation (jl]) represent the Fourier transforms of 
the electronic coupling coefficients and are calculated from the energy in equation ([3]), 
i.e. 

, ^ d'EjRX) ... 

K denotes the EDF (CF and DF in the present model, see equation ([2])) in an elementary 
cell. The B coefficients describe the coupling between the EDF and the displaced 
ions (bare electron-phonon coupling), and the coefficients C determine the interaction 
between the EDF. The phonon frequencies (^^(q) and the corresponding eigenvectors 
e°'{qa) of the modes (qcx) are obtained from the secular equation for the dynamical 
matrix in equation (jlj), i.e. 

J2ttf{^)e^{^)=^'{^)en^)- (7) 

The Eqs. (l4])-([7j) are generally valid and, in particular, are independent of the specific 
model for the decomposition of the perturbed density in equation ([2]) and the pair 
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approximation in equation for the energy. The lenghty details of the calculation of 
the coupling coefficients B and C cannot be reviewed in this paper. They are presented 
in [13]. In this context we remark that the coupling matrix CKK'(q) of the EDF-EDF 
interaction, whose inverse appears in equation (jl]) for the dynamical matrix, can be 
written in matrix notation as 

C = U-^ + V. (8) 

is the inverse of the irreducible (proper) polarization part of the density response 
function (matrix) and contains the kinetic part to the interaction C while V embodies 
the Hartree and exchange-correlation contribution. needed for the dynamical 

matrix and the EPI is closely related to the (linear) density response function (matrix) 
and to the inverse dielectric function (matrix) e~^, respectively. 

Only very few attempts have been made to calculate the phonon dispersion and 
the EPI of the HTSC's using the linear response method in form of density functional 
perturbation theory (DFPT) within LDA [191 [201 [211 [22] . These calculations correspond 
to calculating 11 and V in DFT-LDA and for the metallic state only. On the other 
hand, in our microscopic modeling DFT-LDA-SIC calculations are performed for the 
various densities in equation ([2|) in order to obtain the coupling coefficients B and V. 
Including SIC is particularly important for localized orbitals like Cu3d in the HTSC's. 
Our theoretical results for the phonon dispersion [6l[8l[T6], which compare well with the 
experiments, demonstrate that the approximative calculation of the coupling coefficients 
in our approach is sufficient, even for the localized Cu3d states. Written in matrix 
notation we get for the density response matrix the relation 

C-^ = U{1 + VU)-^ = Ue-\ e=l + VU. (9) 

The CF-CF submatrix of the matrix 11 can approximatively be calculated for the metallic 
(but not for the undoped and underdoped) state of the HTSC's from a TEA of a single 
particle electronic bandstructure. In this case the electronic polarizability H in tight- 
binding representation reads: 

n . 2 >^ /„,(k + q)-/„(k) 

X [C:„(k)an'(k + q)] [C:,MC.'n'{k + q)]* . (10) 

/, E and C in equation ffTOj) are the occupation numbers, the single-particle energies and 
the expansion coefficientes of the Bloch-functions in terms of tight-binding functions. 

The self-consistent change of an EDF at an ion induced by a phonon mode (qa) 
with frequency (^^(q) and eigenvector e°'(qcr) can be derived in the form 



^X-"(q)u,(qa) 



= 5a(q^)e^'*^^ (11) 
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with the displacement of the ions 

h 



1/2 



^ u„(qa)e^'^^^ (12) 

The self-consistent response per unit displacement of the EDF in equation (fTTj) is 
calculated in linear response theory as: 

X(q) = n(q)e-i(q)B(q) = C-i(q)B(q). (13) 

The generalization for the quantity 11 in Eqs. ([8]) and iQ needed for the kinetic part of 
the charge response in the nonadiabatic regime, where dynamical screening effects must 
be considered, can be achieved by adding {hjj+irf) to the differences of the single-particle 
energies in the denominator of the expression for 11 in equation (ITOl) . Other possible 
nonadiabatic contributions to C related to dynamical exchange-correlation effects and 
the phonons themselves are beyond the scope of the present approach. Using equation 
([9]) for the dielectric matrix, e, and the frequency-dependent version of the irreducible 
polarization part, 11, according to equation (fTOll . the free-plasmon dispersion is obtained 
from the condition, 

det[£,«,(q,cu)] =0. (14) 

The coupled-mode frequencies of the phonons and the plasmons must be determined self- 
consistently from the secular equation (J?!) for the dynamical matrix which now contains 
the frequency uj implicitly via 11 in the response function . Such a nonadiabatic 
approach is necessary for a description of the interlayer phonons and the charge-response 
within a small region around the performed in [6l |12] and in the present paper 

for Sr2Ru04. 



3. Results and discussion 



3.1. Modification of the LDA-based handstructure for Sr2Ru04 and Ai phonons 

We shortly recall the construction of a modified LDA-based model for La2Cu04 which 
better describes the real anisotropy of the HTSC's that as a rule is not correctly 
represented in LDA being generally too isotropic. In our effort to obtain a reliable 
description of the Ai phonons polarized and propagating along the c-axis of the cuprates 
[12] we have shown that these modes are highly sensitive with respect to details of the 
c-axis coupling and thus a very accurate representation of the electronic /c^-dispersion is 
needed. Typical LDA based computations lead to an overestimation of the /c^-dispersion 
and give imprecise results for the Ai modes. 

In figure [H we illustrate the degree of inaccuracy for the Ai phonons, in La2Cu04 
as obtained in [T2j. Figure [It^a) shows the result of an LDA-based tight-binding 
representation of the electronic bandstructure (31 band model, 31BM) and figure [T](b) 
gives the outcome of a modified bandstructure (M31BM) with reduced first neighbour 
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Figure 1. Calculated phonon dispersion for La2Cu04 of the c-axis polarized Ai modes 
[121 based on the 3 IBM for the proper polarization part 11^^' (a) and the M31BM (b), 
respectively. The experimental data are represented as open squares (□). The full dot 
(•) denotes the Of mode and the open circle (o) the ylf„(ferro) mode. 



Oj^y-La parameters of the 31BM by 1/6 and of the first neighbour La-La parameters by 
1/3. As seen from figure [1] this leads to a much better result for the phonon modes. 
The characteristic experimental features of the dispersion i.e. the step-like structure of 
the second highest branch and most significant the third highest branch with the steep 
dispersion towards the Z point are not well reflected by the calculation with the 31BM 
as input for the proper polarization part in equation (fTOl) . Using, however, the M31BM 
the calculated phonon dispersion is in good agreement with the experiment as can be 
seen in figure [Tt^b). The characteristic features are now well described. 

We also find as a consequence of the increased anisotropy in the M31BM a 
rearrangement of the three Z point modes with the highest energy. While in the 31BM 
the strongly coupling Of mode is the second highest mode and the y4f^(ferro) mode the 
lowest one in the M31BM Of is the lowest and 742u(ferro) the second highest of the three 
modes. Thus Of is the end point of the steep branch. During this mode rearrangement 
the frequency of Of stays nearly constant while that of Af„ (ferro) is strongly increased. 

We have recalled these facts for La2Cu04 because for Sr2Ru04 a comparable high 
sensibility of the Ai modes in dependence of the electronic fc^-dispersion occurs and 
an analogous mode behaviour for Of and Af^ (ferro) results in response to a weaker 
fc^-dispersion. Likewise as in La2Cu04 this sensitivity of certain c-axis phonons is used 
to construct an accurate tight-binding description of the electronic bandstructure (BS) 
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Figure 2. Calculated phonon dispersion for Sr2Ru04 of the Ai modes based on 
different models for the electronic bandstructm'e as explained in the text. 27BM (a), 
M27BM1 (b) and M27BM2 (c). The experimental values from [TT] are represented as 
open squares (□). The full dot (•) denotes the Of mode and the open circle (o) the 
Af„(ferro) mode. 



in Sr2Ru04. 

As a first approximation for tlie BS we employ a tiglit-binding representation of tlie 
first principles linearized - augmented - plane-wave (LAPW) BS as obtained within the 
framework of DFT-LDA [23]. This analysis leads to a 27 band model (27BM) including 
Sr4d, Ru4d and 02p states. 

The associated calculated phonon dispersion for the Ai modes of Sr2Ru04 along 
the A ~ (0, 0, 1) direction is depicted in figure[2](a) and for the relevant Ai and A4 mode 
along the A ~ (1, 0, 0) direction in figure[3](a) together with the experimental data points 
[11]. Figures [2](b), (c) and figures [3](b), (c), respectively, display our results for two 
modified 27 band models (M27BM1 and M27BM2) with the tight-binding parameters 
important for c-axis is coupling {Oxy-Oz, Ru-Sr, O^^^-Sr, Ru-O^) reduced by 1/2 and 
1/5 with respect to the LDA-based 27BM. This leads as shown in figures |l](a)-(c) to a 
much stronger anisotropy as in the 27BM, i.e. to a strongly reduced fc^-dispersion which 
is also much weaker than in La2Cu04 |12j . 

In figure H] the reduction of the electronic fc^-dispersion is shown along the cut of 
the Fermi surface (FS) highlighted in figure [5](a) when going from the 27BM to the most 
anisotropic M27BM2 via an intermediate c-axis coupling in M27BM1. The best result 
for the phonon modes in figures [2] and [3l which are less well described in the 27BM, are 
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Figure 3. Calculated phonon dispersion for Sr2Ru04 of a Ai mode ( — ) and a A4 

mode ( ) based on the 27BM (a), the M27BM1 (b) and the M27BM2 (c). The 

experimental data points TT] are characterized by open squares (□) for Ai and by full 
squares (■) for A4. The full dot (•) and the open circle (o) at the Z point represent 
of and Af„, respectively. 
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Figure 4. Electronic fc^-dispersion for Sr2Ru04 along the cut of the Fermi surface 
highhghted in figure Efa) by the black bar. 27BM (a), M27BM1 (b) and M27BM2 
(c). The corners of the surfaces are A = (0.2,0.2,0)^, B = (0.45,0.45,0)^, 
C= (0.2^1 0.2^,^), 1)= (0.45^,0.45^,21). 
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obtained for the highly anisotropic M27BM2, i.e. a nearly two-dimensional BS reflected 
by an even considerably weaker /c^-dispersion as in case of La2Cu04, see figure E] and 

It is enlightening to point out that in context with the Ai and A4 branches in 
figure [3] an investigation of |Hj exhibits that the dispersion cannot be described within 
a normal ionic shell model extended by homogeneous electron gas screening to simulate 
the metallic character of Sr2Ru04. Two specially adapted force constants between Ru- 
O2 and Ru-Oxy had to be introduced to mimic the observed dispersion. Of course these 
force constants do not have an intrinsic physical meaning. On the other hand, from our 
calculation we can conclude that on a microscopic level an accurate electronic dispersion 
along the c-axis is essential to understand the observed mode behaviour. 

For a more detailed discussion of the rearrangement of the phonon modes in figures 
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Figure 5. Comparison of the calculated multi-slieet Fermi surface of Sr2Ru04 in tlie 
fcz = plane in model M27BM2 (a) and the 27BM (b) with the measured Fermi surface 
from 24 . The open dots (o) give the experimental results. For the explanation of the 
small black bar in (a) see figure [H 




Figure 6. Electronic bandstructure En{k) of Sr2Ru04 in the M27BM2 taking into 
account the strongly enhanced anisotropy of the real material as compared with the 
27BM as a typical DFT-LDA based model. 
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Figure 7. Displacement patterns of certain phonon modes relevant for the discussion 
of the phonon dynamics in Sr2Ru04. From left to right we display in the first row 
A^„(ferro), A\^{'\i), O^; in the second row ^f„(Ti), Of, Af„(ferro) and in the third 
row Ai/2 and . 

[2] and [3] as a result of a modification of the fc^-dispersion of the BS it is useful to display 
in figure [7] the displacement patterns of some phonon modes relevant for our studies. 

In the 27BM the frequency of Of is nearly degenerate with y4f^(ferro). We denote 
the latter mode as "ferroelectric" because its displacement pattern looks similar to 
that of the ferro-like A^^{ieiio) mode at the F-point, see figure [3 Here the oxygen 
anions vibrate coherently against the cations in the lattice and as a consequence the 
electric dipole moments generated by the motion add constructively. With decreasing 
/c^-dispersion of the BS in M27BM1 and M27BM2 achieved by a reduction of the relevant 
tight-binding parameters mentioned above by 1/2 and 1/5, respectively, the frequency of 
y4f^(ferro) increases by about 2 THz while Of remains virtually unchanged. This highly 
sensible behaviour with regard to c-axis coupling in Sr2Ru04 has also been demonstrated 
for La2Cu04 [12] where the frequency of (ferro) increases by more than 3 THz as 
can bee read off from figure [H 
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Table 1. Magnitudes of the charge fluctuations \6(^k\ in units of 10^"^ particles from 
equation pT|) excited on the Ru4d, Oxy2p, Oz'^.^i and Sr4d orbitals in the Of and 
Afjj(ferro) mode. The results are given for comparison for the M27BM2 and the 
27BM, respectively. 

Of Afjferro) 
Ru4d 0:,y2p 0^2p Sr4d Ru4d 0^y2^ 0^2p Sr4d 

M27BM 26.98 7.62 0.02 0.49 3.22 0.34 

27BM 27.34 6.43 0.63 3.57 9.19 4.69 

The physical origin for the sensitivity derives from the fact that for symmetry 
reasons in case of y4f„(ferro) CF's can be excited in the screening process only at the 
ions in the ionic layers, i.e. O^, La in La2Cu04 and O^, Sr in Sr2Ru04. The strength 
of these CF's on the other hand is governed by the matrix elements of the proper 
polarization part n^K' in equation ffTOj) of the corresponding out-of-plane ions and, thus, 
depends critically on the coupling along the c-axis. Contrarily, the renormalization 
of Of is nearly exclusively given by the CF's on the Ru4d and Oa;y2p orbitals in the 
plane. For quantitative results compare with the CF's listed in table [H This explains 
the inertness of Of and the sensibility of y4f^(ferro) with respect to the change of the 
electronic fc^-dispersion. 

As a further example for the critical interrelation of c-axis phonon dynamics and 
electronic c-axis coupling we note that the steep branch connecting y42„(ferro) with 
of in La2Cu04 is missing in Sr2Ru04. In the latter case due to mode rearrangement 
y42„(ferro) as the lowest of the three highest modes at F is connected to Af^ (||) that 
is the lowest of the three highest modes at Z and not to Of being the second highest 
mode. 

The BS underlying the M27BM2 is illustrated in figure [61 As already mentioned 
the dispersion of the bands along the A (A;^) direction is extremely small, considerably 
smaller than for the 27BM and for the case of La2Cu04. Thus Sr2Ru04 is a nearly two- 
dimensional Fermi liquid. However, as our calculation in figures |2] and [3] have shown the 
remaining weak three-dimensionality is crucial to obtain a reliable description of certain 
c-axis phonons and the c-axis charge response in general. 

Strong hybridization between Ru4d and Oa;j/2p states is evident from the BS in 
figure [6] because of the striking in-plane dispersion of the bands. Moreover, we find a 
similar shape of the partial density of states (PDOS) of the Ru4d and the orbitals 
around the Fermi energy sp- The states sX ep are more strongly of Ru4d type with some 
admixture of Oxy^^)- The contribution of Oz^^) and Sr4d around is very small. In the 
M27BM2 we obtain for the PDOS at ep ^Ru4d {sp) = 3.375 eV'^ , Zo,,2p (sp) = 0.795 
eV-\ Zo^2p{ep) = 0.002 eV'^ and Zsrid{ep) = 0.001 eV'^ 

There are qualitative differences concerning the origin of hybridization comparing 
Sr2Ru04 with La2Cu04. Hybridization in Sr2Ru04 is strongly favoured because the 
radial extent of the 4d wave functions in Ru is much larger than for the 3d wave functions 
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in Cu. Concurrently this leads to a weaker influence of electron correlation effects, e.g. 
a smaller on-site Coulomb repulsion U for Ru4d compared to Cu3d. In our calculations 
we find for [/(Cu3d) = 1.005 dRy and for ?7(Ru4d) = 0.619 dRy. On the other hand, 
hybridization is promoted in La2Cu04 as compared to Sr2Ru04 because the energy levels 
of the copper and oxygen ions are close by in La2Cu04 but far off between ruthenium 
and oxygen in Sr2Ru04. 

A further remarkable difference of the BS in Sr2Ru04 and La2Cu04 can be read off 
from figure[6]and the BS for La2Cu04 in [12]. In Sr2Ru04 up to three bands close to each 
other are crossing the Fermi level giving rise to a three-sheet FS while in La2Cu04 only 
one band is crossing and other bands are more far away. Thus we have a completely 
different situation as interband transitions are concerned which are on a much lower 
energy scale in Sr2Ru04 as in La2Cu04. So we find for -E„'(k + q) — -En(k) with n' ^ n 
for q at the Z point in case of Sr2Ru04 a minimum of 20.67 meV in the M27BM2 
around k = 0.68^(1,1,0), see also figure [6], while in La2Cu04 we obtain for the Z 
point at k = 0.42^(1,1,0) a minimum of 452.71 meV. This, of course, has important 
consequences for the contribution of the interband transitions to II^k' in equation (fTOl) 
concerning magnitude as well as energy scale. 

For example the low energy scale of the interband transitions in Sr2Ru04 can 
generate damping of possible plasmons already at very low energy. Such a damping 
already present in the coUisionless regime together with additional damping due to 
interactions between the electrons very likely leads to overdamping of the very low lying 
plasmons in Sr2Ru04 along the c-axis. The energy of the latter in the coUisionless 
regime is only 5-6 meV along the A direction as will be shown in section 13. 3[ Thus, the 
degree of the remaining weak three-dimensionality quantified by our computations is 
also very important for a possible existence of plasmons and coupled phonon-plasmon 
modes around the c-axis. Such modes have been predicted for La2Cu04 in our recent 
calculations [12] which, however, is not so anisotropic as Sr2Ru04 and the plasmons 
along the c-axis are at higher frequencies that overdamping can be avoided. 

In figure [5](a) and (b), respectively, the calculated FS in the fcz = plane is shown 
for the M27BM2 and the 27BM, together with the experimental data points ^24j. The FS 
consists of the sheets a, /3, 7. The a sheet is holelike and /3 and 7 are electronlike. The 
7 sheet is dominantly derived from the Ru4da;y orbitals, the a and (3 sheets are primarily 
related to the Ad^z and 4dyz orbitals of Ru which exhibit anticrossing behaviour along 
the zone diagonal. The large electronlike 7 sheet passes right through the hybridization 
gap between the a and [3 sheets. 

The experimental results are not well described by the 27BM in figure [5]^b). In 
particular the crossing between the (3 and 7 sheet in the 27BM which is typical forLDA 
calculations is experimentally not reproduced. This is really a crossing and not an 
anticrossing as for /c^ = the bands related to dxy and dxzl(^yz oribitals have different 
symmetry and no mixing is allowed. On the other hand, our modified BS model 
M27BM2 in figure [5](a) is in good agreement with the experimental FS. 

In order to obtain also a more global impression of the magnitude of the enhanced 
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Table 2. Calculated data for the Fermi surface parameters (Drude plasma energy 
tensor, Fermi velocity tensor) and the anisotropy ratio A^^ = ^p,xxl^p,zz] = 
(«F,a;a;)^^V(«F,zJ^^^ in the 27BM and M27BM2 for Sr2Ru04. For comparison the 
corresponding results are also given for La2Cu04 [12]. 

27BM 1045.94 65.50 6.63 0.88 15.97 7.53 

M27BM2 1017.40 3.57 6.03 0.08 284.99 75.38 

M31BM 648.60 25.25 2.99 0.11 25.69 27.18 



anisotropy in the M27BM2 as compared to the 27BM we compare some FS parameters 
being important for transport hke the Drude plasma energy tensor and the Fermi velocity 
tensor. The Drude tensor ist defined as 

hVLp^ij = I ^6{En{k) - EF)Vkn,i'Vlin,j j (15) 

\ ^ kn / 

and the Fermi velocity is given by 

= (lj2^iEnik) - Ep)v^^,v^^A , (16) 



N 

kn 



with 



The output of our calculation for the 27BM and the M27BM2 is given in table [21 
As can be seen we obtain a large enhancement of the anisotropy ratio for the plasma 
frequencies Afi and of the Fermi velocities A^j^, by about a factor 18 and 10, respectively, 
in the M27BM2 as compared to the LDA-based 27BM. 

Additionally, we have included in table [2] the corresponding calculated data for 
La2Cu04 [12] • Again we recognize the by far larger anisotropy of Sr2Ru04 also in the 
transport properties. The values for ^Ip^zz and (vp^^)^^"^ are significantly increased in 
La2Cu04. 



3.2. Phonondynamics in Sr2RuO/^ - Comparison with La/zCuO^ 

In this section we investigate the phonon dynamics of Sr2Ru04 in the main symmetry 
directions A ~ (1,0,0), S ~ (1,1,0) and A ~ (0,0,1) and continue our discussion 
from section 13.11 of a comparison of important characteristica as found in our earlier 
calculations for La2Cu04. 

For a definite investigation of the phonon renormalization induced by the nonlocal 
EPI effects of DF and CF type mediated by the second term in equation @ a 
quantitative reference model for the calculation of the phonon dispersion not including 
the nonlocal screening effects but representing the important ionic component of binding 
in the material is needed. Such a model sketched in section [2] including approximately 
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Figure 8. Calculated phonon dispersion of Sr2Ru04 for the RIM in the main 
symmetry directions A ~ (1,0,0), S ~ (1,1,0) and A ~ (0,0,1). The various 
symbols representing the experimental results are from [11] and indicate different 
irreducible representations. The arrangement of the panels from left to right 
according to the different irreducible representations is as follows: | Ai | A2(- • • , o), 
A4(— ,n)|A3|Si|S2(---,o), E4(— ,n)|S3|Ai(— ,□), A2(---,o)|A3|. 

the local EPI effects is provided by the ab initio RIM extended via covalent ion softening 
and scahng of certain short-ranged pair potentials. 

The result of the phonon dispersion along the main symmetry directions is shown 
for the RIM in figure [8] and compared with the experimental data from inelastic neutron 
scattering (INS) [TT]. The static effective charges found for the model are Ru^'^^, O'^y'^^) 
0~^'^^ and Sr^'^+. In particular the charges of Ru and Oxy differ considerably from their 
nominal values Ru^^ and O^^. This is due to the strong hybridization of the Ru4d and 
Oxy2]) states that reduces the amplitude of the static effective charges in a mixed ionic- 
covalent compound like Sr2Ru04 because of the charge transfer (CT) from the cations 
to the anions is not complete as in the entirely ionic case. Consistent with the result of 
the effective ionic charges is an enhanced covalent character of the Ru-Oj,^ plane with 
a dominant covalent scaling of the Ru-Oa;^ potential and a more ionic character of the 
Sr-Oz layer with an ionic scaling of the Sr-Oz potential as found in the calculations. 

Altogether the RIM with these modifications yields good structural data for 
the energy-minimized configuration. In detail we get for the planar lattice constant 
a = 3.834A, for the lattice constant along the c-axis c = 12.625A and for the internal 
position of the and the Sr ion z{Oz) = 0.163c and z{Sy) = 0.142c. The experimental 
values at 15 K are a = 3.862A, c = 12.723A, ^(O^) = 0.162c, z{Si) = 0.147c [25]. 
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From these data we extract that the tetragonal distortion of the RuOg octahedra is 
significantly smaller than for the CuOg octahedra in La2Cu04. The ratio of the in-plane 
to out-of-plane distance of the oxygen is 0.9196 in Sr2Ru04 but 0.7706 in La2Cu04. 

The phonon dispersion of the RIM is in reasonable agreement with the experiment 
as far as the branches with lower frequencies are concerned. This points to the 
importance of the ionic component of binding in this material. Large deviations are 
observed for the modes with higher frequencies. In particular the high-frequency OBSM 
Ai/2, Og and Of, see figure [3, are not well described. The corresponding frequencies 
are overestimated as compared with the full calculation including DF's and CP's by 
about 5.8 THz, 7.4 THz and 6.7 THz, respectively. This can be attributed according 
to our approach to the missing screening of the Coulomb interaction by DF's and most 
importantly by CP's. 

Noticeable is the soft S4 mode at the X point and the very low frequency 
of the branch in the measurements (the notation for the irreducible representations 
characterizing the symmetry of the modes S4 and S3 is interchanged with the notation 
of [H]). This mode is associated with the rotation of RuOg octahedra around the 
c-axis. As indicated by our calculations the softness of the rotational mode can be 
considered as a precursor of a structural phase transition which in the meanwhile has 
been observed in the Ca2_a;Sra.Ru04 series [11]. From our calculation in the RIM this 
transition can be expected to be driven essentially by the strong component of the 
ionic forces in Sr2Ru04. Contrarily, Sr2Ru04 does not experience the tilt instability 
related to the S3 branch at the X point found in our calculations for La2Cu04. This 
instability correctly indicates the experimentally observed structural phase transition 
from the high-temperature tetragonal (HTT) to the low-temperature orthorhombic 
(LTO) structure. Also this transition is mainly brought about by the long-ranged ionic 
forces. 

Substituting the smaller isovalent Ca ions for the Sr ions induces a misfit between 
the Ca and Ru ions and generates a close connection between the rotation of the RuOg 
octahedra and the electronic and magnetic properties in the Ca2_a;Sra; RUO4 series [TT] . 
Besides the softening of the rotational mode in the range 0.5 < a; < 1.5 around x ~ 0.5 
a similar structural transition occurs due to the softening of the RuOe tilt mode [26] . In 
this context it is interesting that at x ~ 0.5 in the paramagnetic phase of Cai.5Sro.5Ru04 
the 7 Fermi sheet gains a remarkable fc^-dispersion and thus a more three-dimensional 
character than Sr2Ru04 [27] . This issue also emphasizes an interrelation between the 
electronic fc^-dispersion and the RuOg distortions. 

In figure [9] we display the calculated results of the phonon dispersion admitting 
additionally to the RIM anisotropic DF's. This leads to a better agreement of the 
dispersion curves with the inelastic neutron results. Checking these calculations against 
the RIM the width of the spectrum is reduced towards the experiment. The frequencies 
of the transverse and longitudinal optical modes polarized in the plane are decreased 
by the DF's and thereby the largest LO-TO splitting is reduced at the F-point from 8.8 
THz in the RIM to 6.3 THz when dipolar polarization processes are included. At the 
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Figure 9. Same as in figure[5]but allowing additionally dipole fluctuations in Sr2Ru04. 

Table 3. Calculated charge fluctuations 5Ck in Sr2Ru04 in the M27BM2 excited on 
the Ru4d, O^c'i-V, Oj^2p, 0^2y> and Sr4d orbitals and induced by the OBSM Ai/2, Oj^, 
of and the OX, mode. For the displacement patterns compare with figure [T] Units are 
in 10^^ particles, a positive value means an accumulation of electrons. 

Ru4d 0:r2p Oy2p 0^2p Sr4d 

Ai/2 38.252 4.472 -0.002 

Of 50.875 0.119 0.081 

Of 26.984 7.621 7.621 0.022 0.486 

O^ 12.521 -3.714 -3.714 -2.224 -0.319 

same time the frequencies of the modes at F and Z are lowered considerably by the 
DF's in 2;-direction and the large splitting of the A2„(ferro) mode of 10.13 THz in the 
RIM is reduced to 8.85 THz. 

The high-frequency OBSM Ai/2, Of and Of are also decreased due to the screening 
by the DF's but are still overestimated by about 3.9 THz, 5.1 THz and 4 THz as 
compared with the full calculation displayed in figure [in] including additionally CF's. 
Thus, the softening of these modes compared with the RIM is of the order 2-3 THz and 
the remaining effect of the renormalization can be assigned to the CF's on the Ru4d 
and Oa;y2p orbitals. 

The quantitative results for the CF's according to equation ( ITTl) are listed for the 
OBSM in table [3l The CF's for O^ are also given in the table which prove to be 
very different from those of Of despite the same local displacement pattern. As can 
be extracted from figure [7] for O^ and Of the apex oxygens O^ move locally in phase 
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Figure 10. Same as in figure [5] but including additionally dipole fluctuations and 
charge fluctuations on the basis of the M27BM2 in adiabatic approximation. 



against or away from the RuO layers. Because of the weak screening along the c-axis and 
the corresponding large strength of the nonlocal EPI we can expect these vibrations to 
induce CP's in the RuO layers, see table [3] for the calculated results. However, the CP's 
differ qualitatively between Of and O^, respectively, because for the P-point vibrations 
the following sum rule can be derived for the CP's 5C,k [TO] - 

5^5C(Pa) = 0. (18) 

The sum over k in equation flTSl) runs over the CP's in an elementary cell of the 
crystal and a denotes the various eigenmodes at P. Physically equation (fT8|) means 
that local charge neutrality of the elementary cell is maintained in the long-wavelength 
limit (q — > 0), i.e. the mode induced CT is organized within the cell such that from 
the outside the cell looks electrically neutral. Equation flTSj) puts a strong constraint on 
the possible CT in the cell and thus on the screening by CP's. This restriction explains 
the small renormalization of the frequency of only 0.581 THz of O^ in relation to Of 
where the renormalization due to CP's is 4.046 THz. In [11] the pronounced softening 
and the broadening of Of is attributed to a strong coupling between this phonon and 
an interlayer charge transport. Such an interpretation is supported by our calculations 
for of where the CP's can be read off from table [3] and the CT is illustrated in figure fTTl 
in terms of the nonlocal part of the displacement induced rearrangement of the charge 
density according to 

5p„(r,qa) = 5C(q, ^)p.(r - R^)- (19) 
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Figure 11. Contour plot in the RuO plane (a) and perpendicular to the RuO plane 
(b) of the nonlocal part of the displacement-induced charge density redistribution for 
the of mode calculated according to equation (fT9| for the M27BM2. The units are 
in 10~'*e^/a|. The phase of Of is as in figure [T] Full lines ( — ) indicate accumulation 
of electrons in the corresponding region of space and broken lines (—••■—) indicate 
regions of repulsion for the electrons. 

The CF's in equation (fT9l) are obtained from equation ( fTTl) and the form-factors 
from equation ([2]). 

In contrast to the constraint for expressed by equation f|T8|) no such restriction 
must be satistfied in the metaUic phase for Of. While in only an intracell CT is 
allowed summing up to zero, Of generates CF's of the same sign in the cell (table [3]). 
This finally leads to CF's of alternating sign in consecutive RuO layers (cells), i.e. an 
interlayer (intercell) CT is set up, see also figure [TTl The latter provides an effective 
screening mechanism of the Coulomb interaction and explains the strong renormalization 
of Of as compared with O^. 

In the adiabatic approximation used in the calculations so far the interlayer CT is 
instantaneous. It has to be be replaced by a dynamic collective charge transfer in case 
a phono n-plasmon scenario would be realistic as for La2Cu04 in a nonadiabatic region 
nearby the c-axis. The question if such a scenario is also likely for the more anisotropic 
Sr2Ru04 is discussed in section 13.31 
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Figure 12. Calculated dispersion of the highest Ai and Si branch in La2Cu04 [8J (a) 
and Sr2Ru04 (b). In (b) these branches have been obtained for three different values of 
the on-site repulsion U(Ru4d) of the Ru4d orbitals and are characterized by different 
line types, U(Ru4d)=0.619 dRy (—), U(Ru4d)=0.819 dRy (• • •) and U(Ru4d)=1.005 
dRy ( — ). The experimental values are indicated as open squares (□) [11]. 

Table 4. Calculated magnitudes of the charge fluctuations |5Ck| in units of 10"'^ 
particles on the Ru4d orbitals generated by the Ai/2 and the Og mode in dependence 
of the on-site Coulomb repulsion U(Ru4d). v is the frequency of the mode in units 
of THz. The results are listed from left to right taking for U(Ru4d) 0.619 dRy, 0.819 
dRy and 1.005 dRy. The latter value matches the calculated value of U for the Cu3d 
orbital in La2Cu04 and the first value is the calculated ab-initio result for the Ru4d 
orbitals. 

Ai/2 Oj^ 

|5CRu4d| 38.252 26.603 17.871 50.875 33.764 22.048 
V 19.635 20.642 21.381 21.178 22.875 24.031 

An important difference is found when considering the OBSM in the cuprates and 
in Sr2Ru04, respectively. The anomalous softening of the OBSM phononanomalies, 
see figure [TWa). in particular of the half breathing mode Ai/2 which is typical for 
La2Cu04 [8J, and most probably generic for all the cuprate based HTSC's [9], [28] is 
strongly reduced in our calculation as compared with La2Cu04. In the experiments the 
softening is completely absent, see figure [T2](b). Our calculations show that there is 
a dominant microscopic reason for the vanishing of the OBSM anomalies in Sr2Ru04, 
namely the magnitude of the on-site Coulomb repulsion U(Ru4d) of the 4d orbitals of 
ruthenium. The latter is obviously underestimated in our ab-initio calculation based on 
the ionic form-factor p^- Compared with the calculated value for the on-site repulsion 
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Figure 13. Calculated phonon dispersion in the M27BM2 of the three lowest Ei 
branches displaying a complex anticrossing behaviour. The open squares (□) are the 
experimental data [llj . 

U(Cu3d)=1.005 dRy in La2Cu04 which leads to good results for the frequencies of the 
OBSM in this material, the calculated value for ruthenium, U(Ru4d) = 0.619 dRy, seems 
not so reliable because of the more incomplete 4d shell of Ru and the spherical average 
assumed in our computations for the orbital densities. So, we increased the value of the 
Ru4d on-site repulsion in two steps towards the calculated result for the Cu3d orbital 
for comparison. As a consequence the phonon anomalies completely vanish in these 
calculations in agreement with experiment as displayed in figure [TWb). Simultaneously 
by increasing U(Ru4d) the CF's on the Ru4d orbitals that are a measure for the nonlocal 
electron-phonon coupling are strongly decreased as can be seen in table HI 

In |29] and [30] sharply peaked magnetic fluctuations at qo = 0.6^(1,1,0) have 
been observed for Sr2Ru04. These fluctuations have been related to dynamical nesting 
properties between the fiat a and (3 sheets of the FS, compare with figure [5](a). Braden 
et al. [TT] have performed a search for Kohn anomalies in the phonon spectrum related 
to this nesting structure. Promising candidates would be the high frequency Si modes 
which are coupled to the Ru4d CF's. However no signature for a Kohn anomaly is 
present in the experiments as well as in the calculation in figure [THl Another possibility 
discussed in [11] is the longitudinal acoustic Si branch where a small dip is observed 
just around qo, see figure fT3l From our calculations it seems very unlikely that this dip 
is a nesting effect because our computations of the phonon dispersion shows that this 
feature is well described by a complex anticrossing between the three lowest branches 
with Si symmetry, as shown in figure [131 So we do not find in Sr2Ru04 any evidence 
for EPI driven by nesting of the FS. 

3.3. Search for a phonon-plasmon scenario in Sr2RuOi 

We investigate for Sr2Ru04 the possibility of a phonon-plasmon scenario along the c- 
axis and in a small region around this axis by performing calculations of the coupled 
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electron-phonon dynamics in the nonadiabatic regime. This is achieved approximatively 
by allowing the proper polarization part H^k' of the DRF to depend on the frequency 
of the perturbation. Different from static screening in the adiabatic approximation this 
leads to dynamical screening of the Coulomb interaction. A nonadiabatic treatment of 
the charge response and a related phonon-plasmon coupling has been shown recently 
to be a realistic option for La2Cu04 if the real anisotropy of the material is taken into 
account in the theory [T2] . 

Due to the much weaker fc^-dispersion of the electronic band-structure obtained in 
the computation for Sr2Ru04 as compared with La2Cu04 we can expect much lower 
plasmon frequencies along the c-axis. This conjecture can be shown to be true for 
example by calculating at the Z point the free-plasmon frequency according to equation 
( IT4|) and checking the result against La2Cu04, see figure [TH From the zero crossing of 
the real and the imaginary part of det (-^) '^)] "^^ S^t for Sr2Ru04 in the M27BM2 
the very small value of ujp{Z) = 5.55 meV (1.342 THz) while for La2Cu04 we have 
LJp{Z) = 39.33 meV (9.509 THz) [12\ which is by a factor of seven larger. This also 
expresses the much stronger anisotropy of Sr2Ru04 in relation to the cuprates. 

In the analysis of the phonon-plasmon scenario we have assumed the coUisionless 
regime, i.e. the quasiparticles (QP) do not scatter each other and damping is only 
possible via intraband electron-hole decay and low-lying interband transitions. However, 
in the cuprates and also in Sr2Ru04 QP scattering due to electron-electron interaction 
and other degrees of freedom is certainly important. This leads to damping and possibly 
overdamping of c-axis plasmons at sufficiently low frequencies and so the latter may cease 
to be well defined collective excitations of the Fermi liquid. 

In order to make an estimate of the range of frequencies where c-axis plasmons 
should exist we examine the zero crossing of the real part of the well known dielectric 
function of an electron gas 



£00 is the high-frequency dielectric constant, 7 = 7 the scattering rate of the QP's along 
the c-axis. Equating e^u) to zero we have the equation. 



Using for 7 the value given in {31] in the normal state of La2Cu04 and the calculated 
value for the high-frequency dielectric constant along the c-axis ~ 2 [16j we obtain 
cu^ > for Up > 7.21 THz. For smaller plasmon frequencies we have u"^ < 0, i.e. the 
plasmon has a diffusive pole and hence the c-axis plasmon is overdamped in the normal 
state. It should be remarked that larger values for would increase the critical value 
for Up beyond that the plasmon could exist. For example setting = 3 or = 4 
we would have Up > 8.83 THz or 10.19 THz, respectively. Assuming a similar c-axis 
scattering rate of the QP's for Sr2Ru04 we regard cUp ~ 7 THz as a reasonable limiting 
frequency beyond which plasmons most likely are not overdamped. 



e{uj) = e. 



00 
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Figure 14. Real- and imaginary part of det[£KK' (q, 'j-')] at the Z point in arbitrary 
units for Sr2Ru04 in the M27BM2 (a) and for La2Cu04 in the M31BM 12J (b). The 
real part is given by the full line ( — ) and the imaginary part by the dotted line (•••). 
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Figure 15. Calculated nonadiabatic coupled phonon-plasmon dispersion within the 
M27BM2. (a) Ai modes ( — ) in the coUisionless regime, the coupling modes at T 
(^2ti(ferro, na) and at Z{0^ (na)) are shown as black dots (•). (b) A'^ modes ( — ) along 
the A' = {e^,0, ^) direction from Z (e = 0) to Z'{e = 0.02). The black dot (•) at Z 
and Z' are the Of (na) and approximately the Of mode in adiabatic approximation, 
respectively. In both figures the broken line ( — ) denotes the free-plasmon branch and 
the dotted line (• • •) the borderline for damping due to electron- hole decay. 

Thus, from our calculations of the plasmon along the c-axis (A direction) displayed 
in figure [TSl^ a) which are at very low frequencies around 1 THz we conclude that the 
latter are overdamped. On the other hand at larger plasmon frequencies which result 
in case the wavevector q is not strictly parallel to the c-axis, as shown for example in 
figure [TST b). the plasmon can acquire enough energy to resist overdamping and may 
exist collective excitation. 

In figure [TST a) we identify instead of the six Ai branches of the adiabatic 
approximation (compare with figure [2](c)) an additional branch due to phonon-plasmon 
coupling. The free-plasmon has frequencies around 1 THz and is nearly dispersionless 
along the A direction. This is in contrast with the situation in La2Cu04 [12j where the 
free-plasmon has significant higher frequencies and increases from about 9.5 THz at the 
Z point to about 12.8 THz at F. Thus, the plasmon is not overdamped along the c-axis 
in La2Cu04 while the phonon-plasmon scenario calculated in the coUisionless regime for 
Sr2Ru04 as displayed in figure fTSTa) is not realistic. 

Analogous to the situation in La2Cu04 [HI [12] the calculated free-plasmon 
frequencies rapidly increase in proportion to the transverse component of the q vector 
in a small nonadiabatic region of the charge response around the c-axis. We have 
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performed calculations along the A" = (e^, 0, and h.' = C (^^) O5 T") direction in 
this region, see figure [161 for different small e- values to study a possible phonon-plasmon 
scenario in Sr2Ru04. We do not go into the details of these computations here but what 
is important concerning the question of overdamping is that for e > 0.0040 in case of 
the A' direction and for e > 0.0025 in the A" direction the free-plasmon frequency is 
larger than 7 THz, the estimated value for the phonon-plasmon scenario to exist. From 
the calculations we find that the longitudinal A2„(ferro, na) mode at F (figure fTST a)) is 
virtually unchanged at e = 0.0040 and so does exist around 19 THz. The same holds 
true for Of (na) at the Z point being practically the same as of for e = 0.0025 in the 
A" direction. Similar as for La2Cu04 [12] with increasing e A2u(ferro, na) and Of' (na) 
rapidly leave the spectrum to high frequencies. 

In case of Of (na) this mode behaviour can also be extracted from figure fTST b) where 
we have presented our calculated results of the dispersion of the coupled phonon-plasmon 
A'l modes along the A' = [e^, 0, ^) direction from Z = (O, 0, ^) to Z' = [e^, 0, ^) 
for < e < 0.02. The highest mode at Z is the Of(na) mode also seen in figure [TST a) 
as a black dot. The broken line is the dispersion of the free-plasmon branch calculated 
from equation (HM and the dotted line is the borderline for damping due to electron- 
hole decay investigated from j^™^^ (-S'n(k) — -E„(k + q)) for the bands crossing the Fermi 
level. From this figure the range of the region with a nonadiabatic charge response can 
be estimated. It is characterized by the steep branch which ultimately converges to 
the frequency of the adiabatic Of mode (black dot at Z'). Similar as in La2Cu04 this 
nonadiabatic region is very small and can be estimated from the figure at about e = 0.01. 
Unlike the case of La2Cu04 the phonon-plasmon dispersion in figure [TST b) should not 
exist in the full range of frequencies according to our estimate for overdamping below 
about 7 THz. 

Parallel to the investigation of the Of mode in La2Cu04 P,[l2], the line-broadening 
of this mode observed by [11] in Sr2Ru04 can be understood from the calculated phonon- 
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Figure 16. Schematic representation of tlie nonadiabatic region in the {kx, k^) plane 
with the directions A' = (e^,0, ^), A' = ((e^,0,^) and A" = (e^,0,C^). 
C G [0,1]- na: nonadiabatic region; ad: adiabatic region. 
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plasmon scenario. Experimentally there is a limited wave vector resolution in INS for the 
transverse direction perpendicular to the c-axis which for the experiments in La2Cu04 
is on the average e = 0.03 [32]. Thus, the relevant frequency range sampled in the 
measurement for Of is over the steep branch beyond which the mode develops an 
of -like displacement pattern. According to the calculations this sampling occurs for 
mode frequencies larger than 12.5 THz and leads to a corresponding broadening of the 
linewidth for Of observed in the experiment. 

Because the region with a metallic adiabatic charge response outweights by a factor 
of roughly three the nonadiabatic region we can attribute to the Of mode an adiabatic 
frequency of about 15.3 THz. The A'^ branch starting at Of (na) in figure [TMb) has not 
been observed in the experiments so far because a very high q-space solution transverse 
to the c-axis would be needed to resolve the mode dispersion in the small nonadiabatic 
sector. 

In [12] we have pointed out for La2Cu04 the relevance of a strongly coupling 
nonadiabatic Of mode. So it is interesting to see if such a large EPI also exists 
for Of(na) in Sr2Ru04. As a parameter for the strength of the EPI in a certain 
mode (qcr) we have used in the past the orbital averaged changes of the selfconsistent 
crystal potential 5VK(gb"), [Gj [12]. In case of La2Cu04 we obtain in units of meV 
|(^Vcu3d| = 955.98 and |5Vb2p| = 515.98 which should be compared with the considerably 
reduced strength |(5l^u4d| = 277.24 and |(5Vb2p| = 316.64 in Sr2Ru04. 

Altogether, we conclude from our calculations for Sr2Ru04 that nonadiabatic EPI 
via phonon-plasmon coupling is not possible strictly along the c-axis but most likely 
away from this axis in a very small region at higher free-plasmon frequencies. Moreover, 
the strength is significantly smaller as in La2Cu04. 

We have argued above that the longitudinal ferroelectric A2„(ferro, na) mode exists 
in the nonadiabatic sector around the F point at about 19 THz. This nonadiabatic 
result helps to understand why c-axis optical activity seen in the experiments [33] is 
possible despite the fact that Sr2Ru04 is in the metallic phase. In [33] three c-polarized 

phonon structures at 6.04 THz, 10.88 THz and 14.51 THz have been detected 
which agree well with the calculated transverse A^y^ modes in the M27BM2 at 5.16 
THz , 10.42 THz and 14.49 THz. The by far largest oscillator strength is obtained 
for the mode at 10.88 THz and according to our calculation this is the ferroelectric 
mode A2„(ferro). Consistent with the large oscillator strength of this mode in the 
experiment our nonadiabatic result for the longitudinal ferroelectric mode A^^^ (ferro, na) 
predicts a very large splitting of about 9 THz. This splitting dominates the infrared 
response for polarization along the c-axis. Further we find in the calculations two 
other longitudinal A2U modes at 12.379 THz and 6.964 THz which define corresponding 
splittings with the transverse modes above. 

From a theoretical point of view it is important to remark that optical activity in the 
metallic phase cannot be explained using the adiabatic approximation for a calculation 
of the phonon dispersion as is commonly done by applying static DFT for the metal 
[T9| [20l [2n [22] . In such calculations there will be no LO-TO or A2U splittings being a 
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measure for the oscillator strength because the transverse effective charges vanish if the 
metallic phase is treated in adiabatic approximation and consequently the induced dipole 
moments defining the oscillator strength in the dielectric function (matrix) vanish too. 
Thus, there will be no optical activity in a metal by the phonons treated in adiabatic 
approximation which underlies static DFT. 

4. Summary and conclusions 

We have shown that for a reliable description of the c-axis charge response in Sr2Ru04 
and in particular the phonon dynamics along this axis a typical LDA-based model is 
too isotropic and must be modified to account for the much weaker interlayer coupling 
in the real material. Similar as in the cuprates studied earlier the large anisotropy in 
Sr2Ru04 is considerably underestimated in DFT-LDA calculations of the electronic BS. 

Prom our investigations we conclude that an accurate representation of the very 
fainty A;2;-dispersion of the BS nevertheless is essential to understand the c-axis polarized 
Ai modes. While the LDA-based model (27BM) is insufficient a modified much more 
anisotropic model (M27BM2) is well suited. 

We have examined the multi-sheet FS of Sr2Ru04 with the two models and found 
that the experimental results for the FS are not well described within the 27BM. On the 
other hand, the M27BM2 shows a good agreement with the measured FS. This model 
can be considered to represent the real anisotropy in Sr2Ru04 sufficiently well. 

We have also calculated the magnitude of the strongly enhanced anisotropy in 
the M27BM2 compared with the 27BM more globally in terms of some Fermi surface 
parameters important for transport properties. Checking both models against each 
other we find an enhancement of about a factor 18 for the anisotropy ratio of the 
Drude plasma energy tensor and of about a factor of 10 for the Fermi velocity tensor. 
While the calculations demonstrate that Sr2Ru04 is a nearly two-dimensional Fermi 
liquid significantly more anisotropic than La2Cu04 the remaining very weak three- 
dimensionality is crucial to achieve a solid representation of the c-axis charge response 
and of certain c-axis phonons. 

We have calculated in detail the phonon dispersion in Sr2Ru04 along the main 
symmetry directions in the BZ and compared the result with the dispersion of the 
HTSC La2Cu04 being structural isomorphic. A good overall agreement of the calculated 
results is found for Sr2Ru04 with the INS data. Comparing the results of the calculated 
phonon dynamics in the RIM with a model allowing additional DF's and finally the full 
model including DF's and CF's the renormalization of certain modes which are strongly 
coupled via nonlocal EPI effects of DF- and CF-type has been analyzed in detail. In 
this context an important point of the comparison with La2Cu04 is that the anomalous 
softening of the high-frequency oxygen bond-stretching modes, being generic for the 
cuprate based HTSC's, is strongly reduced or completely absent in Sr2Ru04 depending 
primarily on the magnitude of the on-site Coulomb repulsion of the Ru4d orbitals. 

We have also investigated the possibility of a Kohn anomaly in the Ei modes of 
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Sr2Ru04 as discussed in the literature. In our calculations we do not find any evidence 
for such an anomaly driven by nesting of the FS. Instead, the dip in the lowest Ei 
branch seen in the experiments is well explained in our computations by an anticrossing 
effect of the three lowest Si branches. 

Finally, we have examined the question of a possible phonon-plasmon scenario in a 
small region of nonadiabatic charge response around the c-axis which has been shown 
to be a realistic option in La2Cu04. 

Due to the much weaker electronic /c^^-dispersion obtained for Sr2Ru04 the 
calculated free-plasmon frequencies along the c-axis are about a factor of eight smaller 
in the collisionless regime than in La2Cu04. We have argued that damping generated 
by the interactions between the QP's as well as interband transitions which are at much 
lower energy scale than in La2Cu04 leads to an overdamping of the plasmon strictly 
along the c-axis in contrast to La2Cu04. However, a coupled phonon-plasmon scenario 
becomes likely also in Sr2Ru04 at higher free-plasmon frequencies. This occurs in case 
the wavevector is not strictly parallel to the c-axis but has a small transverse component. 
We find that the strength of this nonlocal nonadiabatic couphng is significantly weaker 
than in La2Cu04. 

Ultimately we have explained by our calculations of the nonadiabatic charge 
response the linewidth of the apex oxygen breathing mode at the Z point and why c-axis 
optical activity as seen in the experiments is possible despite the fact that Sr2Ru04 is 
in the metallic phase. 
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